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ABSTRACT 

The observation that long tracts of RNA are 
associated with replicating molecules of mitochon- 
drial DNA (mtDNA) suggests that the mitochondrial 
genome of mammals is copied by an unorthodox 
mechanism. Here we show that these RNA-contain- 
ing species are present in living cells and tissue, 
based on interstrand cross-linking. Using DNA syn- 
thesis in organello, we demonstrate that isolated 
mitochondria incorporate radiolabeled RNA precur- 
sors, as well as DNA precursors, into replicating 
DNA molecules. RNA-containing replication inter- 
mediates are chased into mature mtDNA, to which 
they are thus in precursor-product relationship. 
While a DNA chain terminator rapidly blocks the 
labeling of mitochondrial replication intermediates, 
an RNA chain terminator does not. Furthermore, 
processed L-strand transcripts can be recovered 
from gel-extracted mtDNA replication intermedi- 
ates. Therefore, instead of concurrent DNA and 
RNA synthesis, respectively, on the leading and 
lagging strands, preformed processed RNA is 
incorporated as a provisional lagging strand during 
mtDNA replication. These findings indicate that 
RITOLS is a physiological mechanism of mtDNA 
replication, and that it involves a 'bootlace' mech- 
anism, in which processed transcripts are succes- 
sively hybridized to the lagging-strand template, as 
the replication fork advances. 



INTRODUCTION 

Chromosomal DNA replication in both prokaryotes and 
eukaryotes occurs via coupled leading- and lagging-strand 
DNA synthesis, although many plasmids and viruses use 
alternative replication mechanisms (1). Vertebrate mito- 
chondrial DNA (mtDNA), typically organized as super- 
coiled circular monomers of M6kb (2), is comparable 
with many extrachromosomal elements of prokaryotes 
(1). The different nucleotide composition of the two 
strands of mtDNA has led them to be defined as heavy 
(H) and light (L). In the 1970s, replication of mtDNA was 
proposed to follow the strand-displacement model, based 
on electron microscopy (EM) (3) as well as mapping of 
mtDNA 5' ends (4,5). The model envisioned replication as 
initiating via synthesis of an extended RNA primer 
mapping to a portion of the major non-coding region, 
which was designated as the origin of H-strand synthesis 
(O h ) (2). After synthesis of approximately two- thirds of 
the nascent H-strand, a specific site (O l , the 'origin' of 
light-strand synthesis) is exposed in the parental 
H-strand, at which point synthesis of the L-strand was 
proposed to initiate. In this model, synthesis of each 
strand was envisaged to occur continuously from a 
single initiation site (2), with completion of the nascent 
L-strand beyond O h occurring after separation of the 
daughter circles. During the prolonged delay between ini- 
tiation of first- and second-strand DNA synthesis, the 
lagging strand template was proposed to remain single 
stranded. 

Although this model has many attractive features, our 
own studies, using both neutral 2D agarose gel 
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electrophoresis (2D-AGE) and EM (6-8), found no 
evidence of the predicted partially single-stranded replica- 
tion intermediates, when DNA prepared from highly 
purified mitochondria was analyzed. Instead, all 
replicating molecules detected were substantially duplex. 
However, we also showed that extensive single-stranded 
regions were generated in many replication intermediates 
when they were treated with RNase H, which degrades 
RNA strands when they are hybridized to DNA (8). The 
inference that mtDNA replication intermediates contain 
long stretches of RNA/DNA hybrid was also supported 
by immunopurification using an antibody specific for such 
polynucleotides (6). Extensive regions of RNA/DNA 
hybrid are not predicted by the strand-displacement 
model, but the observation that they were easily converted 
to partially single-stranded species by treatment with 
RNase H, or during nucleic-acid extraction from crude 
mitochondrial fractions, suggested that the molecular 
species from whose analysis the strand-displacement 
model was inferred were, in fact, a preparative artifact. 

The extended RNA tracts found in replication inter- 
mediates were found to comprise uniquely the L-strand 
of mtDNA (7,8), which has recently been corroborated 
by another group (9). The identification of lengthy 
tracts of RNA hybridized to the parental H-strand, led 
us to propose that Ribonucleotides are Incorporated 
Throughout the Lagging Strand, RITOLS (8), although 
until now we were not able to infer the mechanism or 
significance of this process. In addition, we found 
evidence for a distinct population of mitochondrial repli- 
cation intermediates (mtRIs) that were resistant to both 
RNase H and single-strand nuclease, and so were inferred 
to comprise fully double-stranded DNA (dsDNA) (10,11). 
Although these are usually in the minority (8), they can be 
readily detected under specific experimental conditions 
(10,11); however, their significance remains to be clarified. 
One obvious possibility is that two or more mtDNA rep- 
lication mechanisms can operate simultaneously, on dif- 
ferent template molecules. 

The unidirectional nature of replication and the ap- 
proximate location of the origin are shared between the 
RITOLS and strand-displacement models. The crucial dif- 
ference is that while RITOLS explains the presence of ex- 
tensive tracts of RNA/DNA hybrid in replication 
intermediates, the strand-displacement model predicts 
instead extensive regions of single strandedness. 
Although RITOLS intermediates have been characterized 
using a variety of techniques (6,8,10), two possible objec- 
tions can be raised against the RITOLS concept. First, 
mtDNA replication intermediates containing tracts of 
RNA/DNA hybrid have not previously been shown to 
exist in vivo. Conceivably, they could arise as an artifact 
during nucleic acid extraction, by hybridization of a 
displaced H-strand, as predicted to exist by the strand- 
displacement model, with excess L-strand RNA. Second, 
even if molecules containing tracts of RNA/DNA hybrid 
could be shown to be present in vivo, it remains to be 
demonstrated that they are true intermediates in DNA 
replication, and not merely dead-end products of tran- 
scription and replication. If RNA-containing mtRIs are, 
indeed, present in vivo, then the most urgent questions 



concern the subsequent fate of these species and the 
source of the RNA, which could provide a clear indication 
of the mechanism of RITOLS mtDNA replication. 

To address these issues, we implemented two proced- 
ures. First, we subjected intact cells, isolated organelles 
and tissue homogenates to nucleic acid cross-linking, 
and tested whether replication intermediates containing 
RNA/DNA hybrid were present in the cross-linked 
material. Second, using metabolic labeling of isolated 
mitochondria with both DNA and RNA precursors, we 
tested whether the species containing tracts of RNA/DNA 
hybrid were genuine replication intermediates in precur- 
sor-product relationship with fully replicated mtDNA. 
These experiments confirmed that RITOLS intermediates 
are present in vivo, and are processed further to fully 
replicated mtDNA. The source of the RNA incorporated 
into RITOLS intermediates was found to be preformed 
transcripts, rather than extended RNA primers generated 
as the replication fork progresses. On the basis of these 
findings, we propose a new model for mtDNA replication. 
Under this 'bootlace model', preformed L-strand RNA is 
incorporated in vivo at the mtDNA replication fork, via 3' 
to 5' hybridization with the displaced H-strand, and is 
replaced by lagging-strand DNA in a subsequent matur- 
ation process. 

MATERIALS AND METHODS 

Mitochondrial DNA preparation from cells and tissue 

mtDNA from sucrose-gradient purified mitochondria 
was prepared as described previously (12), except that 
isolated mitochondria were treated with Proteinase K at 
4°C for 10-45min before detergent lysis and successive 
phenol and chloroform/isoamylalcohol (24:1) extractions. 
In some cases, whole tissue DNA was extracted dir- 
ectly from homogenized rat liver using phenol and 
chloroform, after being solubilized with 1% sodium 
TV-lauroylsarcosinate and incubated with lOOug/ml 
Proteinase K on ice for 30min. DNA was isolated from 
human cultured cells as previously described (11). Gel ex- 
traction and analysis of mtRIs from Balb C mouse livers 
was performed as described previously (8). 

In oYganello labeling of mitochondrial nucleic acid 

Labeling of mitochondrial nucleic acids in isolated organ- 
elles was performed by a modification of a previously 
described method (13). Mitochondria from rat liver were 
isolated at 4°C using DNase and RNase free reagents. 
Isolated liver tissue from 3- to 5-week-old female 
Sprague-Dawley rats was minced and homogenized in 
4ml/g of tissue in sucrose tris edta (STE)-buffer [320 
mM sucrose, 10 mM Tris-HCl (pH 7.4), ImM EDTA 
and 1 mg/mL essentially fatty acid-free bovine serum 
albumin (BSA)] using a motorized tight-fitting teflon 
pestle. Cellular debris was pelleted by differential centri- 
fugation 1000^ max for 5min. The resulting supernatant 
was centrifuged at 9000g max for 5min to pellet 
mitochondria. Mitochondria were washed once in STE- 
buffer, pelleted and equilibrated in incubation buffer [10 
mM Tris-HCl (pH 8.0), sucrose and glucose 20 mM each, 
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65 mM D-sorbitol, lOOmM KC1, lOmM K 2 HP0 4 , 05 uM 
EDTA, lmg/mL BSA, 1 mM ADP, MgCl 2 , glutamate 
and malate 5mM each]. In organelle labeling was per- 
formed at 37°C with rotation, using 4mg/mL 
mitochondria in incubation buffer supplemented with 
dCTP, dGTP and dTTP (50 uM each) and [oc- 32 P]-dATP 
(Hartmann, 3000Ci/mmol) at 6.6 nM. RNA was labeled 
using 6.6 nM [oc- 32 P]-UTP (Perkin-Elmer, 12000Ci/ 
mmol), ATP, CTP and GTP (50 uM each), and all four 
dNTPs (50|iM each). Mitochondria were recovered by 
centrifugation and the nucleic acids extracted as described 
previously (12). 

Nucleic acid modification, fractionation and detection 

Psoralen cross-linking: 4-5 mg aliquots of sucrose- 
gradient purified rat liver mitochondria, or 15-16mg 
aliquots of rat liver homogenate were treated with 0 and 
5 jiM psoralen or 0 and 25 uM psoralen, respectively (4, 5', 
8-trimethylpsoralen, SIGMA) in 1ml of 20 mM HEPES 
[pH 7.8], 50 mM EDTA. For cross-linking of mitochon- 
drial nucleic acids in human osteosarcoma or embryonic 
kidney cells, the cells were grown to 80% confluence on 
90 mm plates, washed with phosphate buffered saline 
(PBS) and treated with 25 uM psoralen in 1 ml of PBS. 
After incubation on ice for 2min in the dark, samples 
were irradiated with 120 mJ/cm 2 of UV light of wave- 
length 254 nm, for 15min. After the cross-linking step, 
nucleic acids were isolated from rat liver homogenates, 
purified mitochondria or human cells, as described 
above. Restriction digestions were performed on 2-10 jig 
of purified nucleic acids according to the manufacturer's 
(New England Biolabs) instructions. Nuclease treatments 
were performed before or after restriction digestion and 
ethanol precipitation as follows: RNase H (Promega) 1 U, 
15min at 37°C; Ribonuclease T x (Roche) 100 U in 10 mM 
Tris-HCl (pH 7.4), 40 mM NaCl, 1 mM MgCl 2 for 30min 
at 37°C. 

Neutral 2D-AGE was carried by the standard method 
(14). Briefly, first dimension separation was 30 V for 16 h 
in a 0.4% (w/v) agarose gel, at room temperature, 
followed by a second dimension electrophoresis step of 
260 mA for 6h in 1% agarose, at 4°C. After electrophor- 
esis, gels were dried or alkaline transferred to nylon mem- 
branes (GE Osmonics Inc.), UV cross-linked and 
hybridized to either DNA or RNA radiolabeled probes, 
based on (15) for rat mtDNA, (16) for human mtDNA 
and (17) for mouse mtDNA (see supplementary data for 
probe details). 

Riboprobes were generated from amplified mouse 
mtDNA containing a T7 promoter, using 32 P-CTP or 
32 P-UTP (Hartmann) and an in vitro transcription kit 
(Ambion). L-strand specific riboprobes were 1, nt 
14903-15 339 (cyt b); 2, 13 874-14280 (cyt b, c-tRNA Glu 
and RNA5); 3, 12788-13334 (RNA5); 4, 11 546-11 742 
(tRNA His , t RNA Ser(AGY) and tRNA Leu(CUN) ); 5, 5568- 
6044 (COXI); 6, 4951-5326 (tRNA Trp , c-tRNA Ala , 
c-tRNA Asn , O l , c-tRNA Cys and c-tRNA Tyr ); 7, 4502-4920 
(ND2) and 8, 3253-3682 (ND1). Nucleotide numbers were 
based on the revised mouse reference sequence (17), with 
target transcripts indicated in parentheses; c-tRNAs are 



the complementary sequence of the corresponding 
transfer RNA (tRNA). Primers are listed in supplemen- 
tary Data. Hybridization of immobilized mitochondrial 
nucleic acid to probe was in 0.25 M sodium phosphate 
(pH 7.4), 7% SDS, 10 mM EDTA overnight, at 60°C for 
riboprobes, or 65°C for DNA. Post-hybridization washes 
were lx SSC three times, followed by lx SSC, 0.1% SDS 
twice, each for 20min at 65°C. Filters were exposed to 
phosphorscreens and scanned using a Typhoon™ 
phosphorimager (GE Healthcare). 

Quantitation of relative signals from Southern blotting 
and in organello labeling 

After in organello labeling, restriction digestion and frac- 
tionation by 2D-AGE, mtDNA fragments were 
transferred to nylon membranes (GE Osmonics Inc.), 
UV cross-linked and exposed to phosphorscreens for 
7-14 days. When the radioactive signal had decayed to 
the point where it was no longer discernible above the 
background signal, the membranes were hybridized to 
DNA radiolabeled probes and exposed for 1-2 days to 
determine the level of mtDNA species in the steady 
state. Signal was quantitated using a Typhoon™ 
phosphorimager (GE Healthcare). Equivalent portions 
of the arcs and spots on the 2D gel blots were selected 
for quantitation, always avoiding overlapping areas. 
These and other radioactive signals from in organello 
labeling were expressed as a percentage of the total 
signal from the panel and standardized to that of the 
steady state detected by Southern hybridization. 

Quantitation of mtDNA and RNA levels 

Transcript and mtDNA levels were determined by quan- 
titative PCR, using GFP as external reference (18). 

RESULTS 

Nucleic acid cross-linking confirms the presence of RNA/ 
DNA hybrid in mtRIs 

Psoralen, a photo-activatable cross-linking reagent that 
introduces predominantly interstrand cross-links into 
nucleic acids, has been used previously to study RNA/ 
DNA hybrids created by in vitro transcription of single- 
stranded bacteriophage fd DNA (19). Crucially, cross- 
linking of RNA/DNA hybrids renders them resistant to 
RNase H (an enzyme that specifically degrades the RNA 
component of RNA/DNA hybrids), owing to the 
accompanying change in conformation from B to 
Z-duplex (20,21), as demonstrated for a model substrate 
(Supplementary Figure SI A). Therefore, interstrand 
cross-linking of mtDNA was performed to distinguish 
whether RNA associates with replication intermediates 
in mitochondria in vivo, or, alternatively, whether this as- 
sociation could be an extraction artifact. The rationale for 
this experiment is as follows. If RNA is in hybrid with 
DNA in bona fide replication intermediates in vivo, 
cross-linking should protect it from digestion by RNase 
H. Conversely, if the RITOLS-type intermediates arise as 
an artifact during extraction, due to hybridization of 
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mitochondrial RNA to a displaced parental DNA strand 
after protease treatment, the RNA will remain RNase H 
sensitive after in vivo cross-linking. 

When intact rat liver mitochondria (Figure 1, 
Supplementary Figure SIB) or tissue homogenates 
(Supplementary Figure SIC) were treated with psoralen 
before mtDNA extraction, mtRIs proved markedly more 
resistant to RNase H than those of controls, when 
analyzed by 2D-AGE after digestion with BspHI. Cross- 
linking applied to intact human cultured cells also 
conferred RNase H resistance on mtRIs, and markedly 
improved the preservation of fragments containing the 
origin(s) of replication, which form so-called bubble arcs 
(Supplementary Figure SID). The fact that the RNA 
associated with mtRIs can be cross-linked in isolated 
mitochondria, fresh tissue homogenates and intact cells, 
without the prior removal of protein, strongly indicates 
that the RNA is, indeed, already hybridized to replicating 
molecules of mtDNA in vivo. 

Radiolabeled dNTPs are incorporated into 
RNA-containing replication intermediates in organello 

Establishing the presence of RNA hybridized in vivo to 
mtDNA does not of itself prove that the RITOLS 
species are true intermediates of mtDNA replication. To 
address this issue, we combined 2D-AGE analysis with a 
protocol for in organello labeling of newly synthesized 
DNA in isolated mitochondria (13). In this procedure, 
the mitochondria are not synchronized for DNA replica- 
tion, and so molecules at all stages in the replication cycle 
should become labeled by the supply of a radioactive 
deoxynucleotide precursor. Both nascent DNA strands 
should become labeled, regardless of whether DNA syn- 
thesis is strand synchronous. This should include mol- 
ecules forming a bubble arc, in an origin-containing 
restriction fragment, as well as molecules forming replica- 
tion fork (Y) arcs or slow-moving Y-like (SMY) arcs 
derived from all restriction fragments of the genome. 
Importantly, if RITOLS is a physiological mechanism of 
mtDNA replication, many of the intermediates labeled in 
organello should be subject to modification by RNase H, 
unless psoralen-UV cross-linked before extraction. 
Furthermore, if RNA-containing species are generated 
at the replication fork, they should be detectable by the 
shortest incubation times required to incorporate label 
into replication intermediates. 

To test these predictions, isolated rat liver mitochondria 
were incubated in the presence of radiolabeled dATP, 
after which the extracted nucleic acids were digested 
with BspHI and analyzed by 2D- AGE (Figure 2A). The 
'In' linear fragments of non-replicating DNA seen prom- 
inently in the steady state (Southern hybridizations) 
(species i, Figure 2B) accounted for only a minor 
fraction of the material labeled in organello after a 5-min 
incubation (species i, Figure 2A). In most other respects, 
the major mtDNA species labeled in organello were similar 
to those detected by Southern hybridization (Figure 2A 
compared to B); i.e. [oc 32 P]-dATP was incorporated chiefly 
into those species previously assigned as mtRIs (7,8,23). 
Next, the restriction enzyme BlpI was tested because it 



cuts rat mtDNA into two fragments of 5.6 and 10.7 kb, 
and the short fragment contains the putative origin region 
(O h ), thus allowing the expected bubble arc to be well 
resolved, in a region of the gel distinct from the 
expected SMY arcs. Here again the majority of putative 
RITOLS intermediates that became labeled in organello 
had the same mobility as those detected in the steady 
state by Southern hybridization (Supplementary Figure 
S2A compared to D), and they displayed a pronounced 
sensitivity to RNase H (Figure 2C versus D). Specifically, 
RNase H digestion collapsed the bubble arc and grossly 
modified the SMY arc, similar changes as those observed 
in parallel (or previously) by Southern hybridization 
[Supplementary Figure S2B and E; and (7,8)]. The 
labeled replication intermediates of the bubble and SMY 
arcs were also largely resistant to single-strand specific 
RNase Ti (Figure 2E), consistent with the interpretation 
that much of their RNA content is in the form of RNA/ 
DNA hybrid. Psoralen/UV cross-linking following 
labeling, but before nucleic acid extraction from 
mitochondria, conferred resistance to RNase H 
(compare Supplementary Figure S2A and B to G and 
H), providing further confirmation that the replication 
intermediates detected by direct labeling are essentially 
the same as those present in the steady state in vivo. The 
amount of DNA label incorporated into mtRIs increased 
linearly with time in both rat and mouse liver 
mitochondria, prepared in the same way (Supplementary 
Figure S3A-C). Labeling was 2- to 4-fold more efficient 
with [oc 32 P]-dATP than with [oc 32 P]-dCTP (Supplementary 
Figure S3D), and pre-incubating mitochondria with the 
DNA chain terminator 2', 3'- dideoxyadenosine-5'- 
triphiosphate (ddATP) inhibited the labeling of mtDNA 
(Supplementary Figure S3E and F). 

One other feature of the in organello-labded material 
that distinguishes it from that detected by Southern hy- 
bridization was the near absence of the Y arcs, represent- 
ing fully duplex DNA intermediates (i.e. fully duplex 
forms resistant to RNase H). Such arcs are seen routinely 
in Southern blots of 2D gels of material from solid tissues 
and cultured cells, and have been proposed to be products 
of conventional strand-coupled DNA synthesis (23). 
Standard Y arcs were detected by Southern hybridization 
in material from the isolated mitochondria used for in 
organello labeling (Supplementary Figure 2D-F), 
indicating that their integrity was not compromised 
during the incubations, although they did not label effi- 
ciently in organello (Figure 2A and C); nevertheless the 
standard Y arc was visible after 30 min of labeling 
(Supplementary Figure S3B). 

RNA-containing replication intermediates labeled in 
organello are in precursor-product relationship with 
mature mtDNA 

If the RNA-containing molecules pulse-labeled in 
organello are genuine intermediates in DNA replication, 
they should become converted to fully replicated mtDNA 
during a subsequent chase. To test this, we labeled rat liver 
mtRIs in organello using [oc 32 P]-dATP for 5 min, condi- 
tions sufficient to incorporate substantial amounts of label 
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BspHI digested DNA from gradient-purified rat liver mitochondria 
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Figure 1. Psoralen cross-linking in vivo renders mtRIs resistant to RNase H. (A-D) DNA from cross-linked or control mitochondria were isolated, 
digested with BspHI and then, where indicated, digested further with RNase H, before separation by 2D- AGE. After transfer, the membrane was 
probed with probe rl (see Supplementary Data for probe details). SMYs 1, 2 and 3 are SMY replication fork arcs, or their RNase H-modified 
counterparts. (E) schematic BspHI restriction map of rat mtDNA. To the right of the map are illustrated interpretations of species detected by probe 
rl, including typical intermediates of the three principle SMY arcs produced by BspHI digestion of rat mtDNA. In the case of SMY1 (inset), the 
effect of RNase H is interpreted according to whether the mtDNA was subjected to psoralen cross-linking in the organelle. Black lines represent 
DNA, red lines are RNA, red crosses indicate RNA/DNA psoralen-mediated interstrand cross-links, DNA/DNA cross-links are not shown for 
simplicity; a black cross is indicative of a blocked restriction site on one branch of the Rl, owing to RNA/DNA hybrid. To the right of the gel 
images are merged false-colour images of A and C, and B and D, providing a direct comparison of the effect of RNase H on control and cross-linked 
samples: red — untreated (A and B), green — treated with RNase H (C and D). Species resistant to RNase H appear yellow; those modified or 
generated by RNase H appear red and green, respectively. BspHI was chosen because it produces a series of well-resolved SMY arcs covering almost 
the entire mitochondrial genome. 
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Figure 2. The majority of in organello synthesized mtRIs have the mobility and RNase sensitivity of those present in vivo. After in organello labeling 
of rat liver mitochondria (6.6 nM [a 32 P]-dATP, 5min, 1 mg protein per reaction), mtDNA was extracted, digested with BspHI and subjected to ID- 
AGE. Phosphorimager analysis revealed a series of SMY replication fork arcs (A), substantially the same as those detected by Southern hybridization 
after probing with rl (B). Note that probe rl detects a specific fragment of mtDNA (i), whereas incubation of mitochondria with [a 32 P]-dATP labels 
all fragments of the mtDNA during the reaction (i, ii, hi). (C-E) In organello labeled rat liver mtDNA was digested with Blpl, separated by 2D- AGE, 
transferred to filter-membrane and phosphorimaged. Where indicated, samples were treated additionally, after restriction digestion, with RNase H or 
RNase TV Interpretations of the arcs and spots appear immediately below the 2D gel images. SMY — SMY replication fork arc; b — origin-containing 
bubble arc of the RITOLS type [see (8), y — standard replication (y) fork arcs, e — eyebrow arc (8), mSMY — modified SMY, mb — modified bubble 
(22), me — modified eyebrow, 16kb — late maturation intermediate where RNA remains at the restriction site at nucleotide 513, but not at 11234, 
representative examples are illustrated. Also illustrated are the In linear restriction fragments of 5.6 and 1 1 kb. Not illustrated are nicked (uncut) 
circles (nc). 
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Figure 3. Pulse-chase experiments establish a precursor-product relationship between RNA-containing mtRIs and mature mtDNA. (A, B and D) 
Five-minute pulse-labeling of rat liver mitochondria with 6.6 nM [a 32 P]-dATP, (A, B) followed by chases of 10, 25, 40, 55, 115 or 175 min, with 
660 fiM unlabeled dATP or (D) following pre-incubation without [a 32 P]-dATP for 115 or 175 min. Extracted DNA was BlpI digested, separated by 
ID (A) or 2D- AGE (B, D) and phosphorimaged after membrane transfer. Brackets demarcate the loss of signal on the replication bubble arc. 
(C) Quantitation of label in mtRIs, during chase, based on phosphorimaging. Values on vertical axes are based on signal from in organelle* labeling 
(nascent) normalized to Southern blot (steady state) signal for the different species analyzed. The resulting ratios for the chase samples were 
expressed as relative to no chase (0' chase), set as 1. 



(Figure 2 A, Supplementary Figure S3), then subjected the 
mitochondria to cold chases of 25-1 75 min, using excess 
non-radioactive dATP. ID- (Figure 3 A) and 2D-AGE 
(Figure 3B) were used to identify the species into which 
label accumulated subsequently. Non-linear (fork-con- 
taining) species larger than the restriction fragments of 
mature mtDNA were labeled preferentially in a 5-min 
pulse, but over the course of a 3-h chase most of the 
label accumulated into the linear BlpI fragments of 5.6 
and 10.7 kb (Figure 3A). 2D- AGE analysis confirmed 
that all the mtRIs diminished during the course of the 
chase (the bubble arc, the SMY arc and the final matur- 
ation intermediate (revealed in the digest as a 16.3 kb 



linear species), whereas label increased in linear fragments 
of 5.6 and 10.7 kb (Figure 3B), representing fully matured 
mtDNA. In the case of the arc of replication bubbles, the 
loss of signal advanced from its base to the tip of the arc 
(Figure 3B), consistent with new rounds of initiation 
incorporating chiefly cold dATP, during the chase, and 
arguing against the incorporation being the result of 
DNA repair. Quantitation of the different species 
indicated that the increase in signal for the unit length 
fragments was essentially linear (Figure 3C), implying 
that replication proceeded at a constant rate throughout 
the experiments. The decrease in signal of the mtRIs 
followed a logical pattern, with (origin-containing) 
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bubble arcs being the first to decline, followed by SMY 
arcs, with the final maturation intermediate (16 kb linear 
species still containing RNA at one BlpI site) transiently 
increasing > 1.5-fold during the chase, before declining to 
low levels (Figure 3C). Replication in the mitochondria 
remained active for the duration of the reactions, as 
mtRIs were efficiently labeled by a 5-min pulse of [oc 32 P]- 
dATP after pre-incubations of 2 or 3 h (Figure 3D). 
Intact, unlabeled mtRIs were also preserved over such 
times, based on Southern hybridization (Supplementary 
Figure S4). 

Labeling of mtRIs with an RNA precursor 

Because the RITOLS model of replication predicts that 
the lagging strand template is hybridized to RNA (8), 
ribonucleotide incorporation into mtRIs was tested, in 
organello, by using [oc 32 P]-UTP in place of [oc 32 P]-dATP. 
Although, there was no detectable labeling of mtRIs after 
10- or 30-min incubations with [oc 32 P]-UTP (Figure 4A 
and B), a faint pattern, similar to that seen with [oc 32 P]- 
dATP, was discernible after 1 h of labeling (Figure 4C), 
and by 2 h of incubation with [a 3 P]-UTP, the 



characteristic pattern of mtRIs was evident (Figure 4D). 
During their synthesis, mitochondrial transcripts are 
expected to be sensitive to RNases, and yet the UTP- 
labeled mtRIs were largely resistant to RNase Ti 
(Figure 4E). In contrast, almost all of the mtRIs that 
became labeled with [oc 32 P]-UTP in organello were 
degraded by RNase H (Figure 4F). Because RNase H- 
treated mtRIs labeled with [oc 32 P]-dATP were modified, 
but still detectable (Figure 2D), and the same was true 
of unlabeled mtRIs in the [oc 2 P]-UTP labeling experi- 
ments, as detected by Southern hybridization (Figure 
4F*), the ribonucleotide precursor cannot have been con- 
verted to DNA once inside mitochondria. As a control, 
the efficient synthesis of bulk mitochondrial RNA was 
verified, during 20 min of [oc 32 P]-UTP labeling in organello 
(Supplementary Figure S5A). 

An RNA chain terminator inhibits the incorporation of 
labeled RNA into mtRIs, without blocking mtDNA 
replication 

Two possible models for RNA incorporation during 
mtDNA replication have been proposed (8): one is 



Nucleic Acids Research, 2013, Vol. 41, No. 11 5845 



In organelle* labeling of rat mtDNA 



[ 32 P]-UTP labeling 



Control 



4 



3' dATP 



B/pl 2D-AGE 



B [ 32 P]-dATP labeling 

B/pl 2D-AGE 



55' 



115' 





• 


B/pl 


i 






1 








Control 


3' dATP 



S c 
2 ° 
S-o 

O CL 

.E < 

< y 



1.2 
0.8 
0.4 
0.0 




-5.6 kb linear -O- 1 0.7 kb linear 



25 55 115 175 

Incubation time (min) 



D 1.2 

< _ 

z o 

3 § 08 
0.4 

O 

■- CO 



0.0 



-mRNA — 0— rRNA 



30' 



1h 2h 

Incubation time 



3h 



Figure 5. An RNA chain terminator rapidly inhibits nascent RNA formation in isolated mitochondria but affects mtDNA replication much more 
slowly. (A) 2D- AGE analysis of rat mtDNA labeled for 120 min with [a 32 P]-dUTP after pre-incubation without (control) or with 20fiM cordycepin 
triphosphate for 20 min. The intactness of the mtRIs is confirmed in Supplementary Figure S5C. (B) 2D-AGE analysis of rat mtDNA labeled for 
5 min with [a 32 P]-dATP after pre-incubation without (control) or with 20jiM cordycepin triphosphate for 55 or 115 min. (C) Quantitation of the 
signals from the linear BlpI fragments of mtDNA (panel B and equivalent samples), during 3' dATP treatment. Note the slow but gradual decrease 
in [a 32 P]-dATP incorporation into mtDNA. [a 32 P]-dATP incorporation was standardized to the level of steady-state mtRIs, based on Southern 
hybridization. (D) Steady-state levels, relative to zero time point, of mitochondrial rRNA (12S and 16S) and mRNA (COX2, ND1, ND3, ND5 and 
CYTB), after the indicated times of incubation with 20 fiM cordycepin triphosphate. Measurements were based on Q-PCR values normalized to those 
for mtDNA. 



based on the successive hybridization of mature tran- 
scripts to the lagging-strand template (the bootlace 
model); the other depends on a primase that generates 
long primers. A third possibility would be that ongoing 
transcription supplies the RNA that hybridize to 
replicating DNA molecules as the fork proceeds, yet this 
could only avoid substantial single-stranded regions 
forming if multiple active transcription complexes were 
arrayed on replicating mtDNAs ahead of the fork 
(Supplementary Figure S6). To discriminate between 
these different possibilities, we added the RNA chain 
terminator cordycepin triphosphate (S'-deoxyadenosine- 
S'-triphosphate or 3' dATP) to in organello labeling reac- 
tions. Although the RNA chain terminator can also block 
synthesis of the primers needed for the initiation of 
mtDNA synthesis, the probability of it doing so will be 
much lower than for termination of the RNAs required 
for expression of mtDNA, which are synthesized as near 
genome-length polycistronic transcripts, or the synthesis 
of the RNA equivalent of Okazaki fragments. Hence, 
there should be a large dynamic range in which an RNA 
chain terminator would effectively inhibit synthesis of the 
long tracts of RNA for the lagging strand (whatever their 



unit size), while having little impact on primer synthesis at 
the origin. Incubation of mitochondria with cordycepin 
triphosphate, before the addition of [oc 32 P]-UTP, 
produced a marked inhibition of free transcripts labeled 
in organello at a drug concentration of 20|iM 
(Supplementary Figure S5A). Cordycepin triphosphate 
also inhibited [oc 32 P]-UTP labeling of mtRIs (Figure 5A 
Supplementary Figure S4B), yet the effect of the inhibitor 
on DNA labeling (with [oc 32 P]-dATP) was much slower 
and less marked (Figure 5C, compare with Figures 2C 
and 5B) and produced only minor changes in the 
relative labeling of different mtRIs. Note that, in 
contrast to its marked effect on ongoing transcription 
(Supplementary Figure S5A), incubation with cordycepin 
triphosphate over 2h did not substantially affect the 
steady state level of mitochondrial RNAs (Figure 5D). 

These results strongly suggest that preformed tran- 
scripts are the source of the RNA incorporated into rep- 
lication intermediates, rather than de novo synthesis via a 
primase, or ongoing transcription. To analyze directly the 
nature of the RNA incorporated into replicating mol- 
ecules, we eluted mouse mtRIs found in the SMY region 
of 2D gels, and analyzed the structure of the RNA therein. 
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We chose the SMY region because, in the Bell digest, it is 
well separated on 2D gels from all other material. As a 
control, we first confirmed that there was no detectable 
trace of even the largest processed mitochondrial RNA 
species resolving in the vicinity of the SMY arcs when 
purified RNA was independently fractionated by 2D- 
AGE (Supplementary Figure S7). The RNA and DNA 
species found in the SMY arc-derived mtRIs were then 
analyzed by blot hybridization, following different enzym- 
atic treatments, heat denaturation and fractionation by 
ID-AGE (Figure 6, Supplementary Figure S8). A series 
of L-strand specific riboprobes revealed the SMY-derived 
mtRIs to contain mainly discrete RNA species, the most 
abundant of which correspond in size to previously 
characterized fully processed transcripts (24-27) from 
around the mitochondrial genome (Figure 6, 
Supplementary Figure S8). Thus, these processed mito- 
chondrial transcripts are proposed to be the source of 
the lagging strand RNA that underpins RITOLS replica- 
tion. It is noteworthy that they included antisense tRNAs 
(Supplementary Figure S8) that have no other predicted 
function than to cover the lagging-strand template, during 
the interval between the initiation of first- and second- 
strand mtDNA synthesis. 



DISCUSSION 

The application of psoralen/UV cross-linking and in 
organello labeling provide compelling evidence that 
RNA/DNA hybrids form naturally during the process of 
mammalian mtDNA replication, and are then replaced 
with DNA, via a maturation process. Furthermore, 
combined with the analysis of RNA species found in 
mtDNA replication intermediates by strand- specific 
probing, they support an unorthodox mechanism for the 
incorporation of RNA into replication intermediates, via 
y to 5' hybridization of preformed L-strand RNAs with 
the parental H-strand, as the replication fork progresses. 
This 'bootlace mechanism', illustrated in Figure 7, 
proposes the successive threading of processed transcripts 
onto the lagging-strand template, to create RNA/DNA 
hybrid tracts in which the RNA strand is discontinuous. 
At most points during replication there will be a portion 
of the RNA that has yet to be hybridized (Figure 7), and 
this might explain why RNase Ti enhanced the signal of 
the SMY arc (Figure 2E). 

Previous EM data showed that the mtDNA replication 
intermediates detectable in steady state are essentially 
duplex over their entire length (6). Furthermore, 
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2D-AGE, enzymatic digestion and immunopurification 
showed that a prevalent class of these duplex intermedi- 
ates contained tracts of RNA/DNA hybrid, implying the 
incorporation of lagging-strand RNA across the whole 
genome, RITOLS (6-8). While we could not exclude 
that these molecules might contain one or more short 
singe-stranded gaps, their overall structure was fundamen- 
tally inconsistent with the predictions of the strand-dis- 
placement model. Moreover, species predicted by the 
strand-displacement model were completely undetectable 
in mtDNA prepared from highly purified mitochondria. 
Nevertheless, these earlier studies left open the possibility 
that the RNA-containing intermediates arose from 
products of strand-displacement replication hybridizing 
with excess mitochondrial RNA during extraction. They 
also did not exclude the possibility that the RNA-contain- 
ing species might be dead-end products. The data pre- 
sented here refute these ideas. Cross-linking in vivo 
(Figure 1, Supplementary Figure SI) protected the 
RNA/DNA hybrid-containing species from degradation 
by RNase H, indicating that they are already present 
before extraction. Metabolic labeling in organello con- 
firmed that they are genuine intermediates in DNA repli- 
cation, as they are in precursor-product relation to mature 
mtDNA (Figures 2 and 3). The RNA found hybridized to 
DNA in replication intermediates corresponds to the 
L-strand of mtDNA (Figure 6, Supplementary Figure 
S8). Hence, we conclude that RNA is laid down on the 



lagging strand as leading-strand DNA synthesis 
progresses around the circular genome. The predicted 
intermediates of the strand-displacement model were not 
synthesized in organello, even under the shortest labeling 
times required for detectable amounts of radiolabel to be 
incorporated into mtDNA. Thus, we can exclude the pos- 
sibility that replication proceeds initially via the strand- 
displacement model, with RNA incorporated into replica- 
tion intermediates in a subsequent step. 

Comparison of the effects of DNA and RNA chain 
terminators (Figure 5, Supplementary Figures S3 and 
S5) showed that concomitant RNA synthesis is not 
required for RITOLS replication. The fact that over 1 h 
of exposure to cordycepin triphosphate failed substan- 
tially to modify mtDNA replication (Figure 5B), while 
effectively inhibiting transcription within 10 min, refutes 
the idea that the lengthy tracts of RNA incorporated on 
replicating molecules of mtDNA are the products of 
ongoing transcription. Equally, the lack of any effect of 
the RNA chain terminator on mtDNA replication 
excludes the possibility of a primase being responsible 
for the synthesis of the lagging-strand RNA. If such a 
primase were sensitive to inhibition by cordycepin triphos- 
phate, replication would arrest, or the lagging-strand 
template would be left bare, profoundly altering the struc- 
ture, mobility and enzymatic sensitivity of the mtRIs [see 
(6,8) for details]. Conversely, a primase insensitive to 
cordycepin triphosphate (28) should incorporate 
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[a P]-UTP into mtRIs irrespective of the presence of the 
RNA chain terminator, which was not the case (Figure 
5A, Supplementary Figure S5B). Thus, the RNA species 
found in RITOLS replication intermediates must logically 
be synthesized before the onset of DNA replication. The 
obvious source of such RNAs is processed (mature) tran- 
scripts, which was confirmed by our analysis of the species 
found in mtRIs isolated from the SMY-arc region of 2D 
gels (Figure 6, Supplementary Figure S8). Although these 
RNA species display some heterogeneity, they are of sub- 
stantial length, and the largest and most abundant corres- 
pond in size with previously characterized mature 
transcripts from the various regions of the genome 
interrogated. This finding provides strong support for a 
modified version of the bootlace model (8), as illustrated 
in Figure 7, in which processed lagging-strand transcripts 
are successively incorporated into replication intermedi- 
ates as the replication fork advances, and are subsequently 
fragmented or recovered during the maturation process by 
an as yet unknown mechanism. Accordingly, the matur- 
ation process might account for some of the aforemen- 
tioned length heterogeneity of the RNAs associated with 
mtRIs. 

One curious observation is that while replication inter- 
mediates composed fully of dsDNA are readily detectable 
by Southern hybridization (10,11,23,29), they were barely 
detectable in the material labeled in organello using DNA 
precursors (Figure 2). As previously suggested, they might 
reflect a parallel mechanism of conventional strand- 
coupled replication occurring alongside RITOLS in a 
separate population of molecules. The fact that they do 
not become labeled in organello, as readily as RITOLS 
intermediates, implies that conventional strand-coupled 
replication is slow in mitochondria, or the isolated organ- 
elles are unable to sustain this mode of replication at the 
in vivo level, perhaps because of the disruption of the mito- 
chondrial network. Alternatively, the intermediates 
composed fully of dsDNA could represent molecules in 
which replication has stalled, and in which, although the 
RNA has been replaced by DNA, their further processing 
occurs only slowly, causing them to accumulate over a 
period of many hours, which is barely detected in 
labeling reactions of maximally 3h (Supplementary 
Figure S3B and C). 

Many new questions arise from this work. The identity 
of the various enzymes and factors that execute and 
regulate the process of RNA incorporation during 
mtDNA replication, as well as the pathway by which the 
RNA containing intermediates are processed subsequently 
during maturation to dsDNA, remain to be elucidated, 
including the extent to which the incorporated RNA is 
used for priming. Poly(A) tails present on processed tran- 
scripts would remain unhybridized to the DNA template, 
but should not constitute a problem for incorporation, 
and may even facilitate the delivery and threading of suc- 
cessive RNA segments as the replication fork progresses, 
or may be used to initiate the maturation process. RNA 
incorporation is likely to require one or more RNA 
helicases to unwind the extensive secondary structure of 
RNAs, especially ribosomal and transfer RNAs. Extensive 
modifications to tRNAs might impair incorporation, and 



so it will be interesting to determine whether RITOLS 
intermediates contain fully or partially unmodified 
tRNAs. 

The evolutionary rationale for the emergence of 
RITOLS, operating via the bootlace mechanism, is not 
obvious. However, for coating of the displaced parental 
strand during DNA synthesis, RNA has a major theoret- 
ical advantage over the mitochondrial single-stranded 
DNA-binding protein, which is the alternative according 
to the strand-displacement model. Although the fidelity of 
RNA synthesis is lower than DNA synthesis, RNA and 
DNA contain the same genetic information. Hybridized 
RNA thus provides an informational back up that could 
be used to repair damage that occurs to the lagging-strand 
template during mtDNA replication, in an environment 
rich in DNA-damaging oxygen radicals. In fact, RNA in- 
corporation affords the possibility of having even less 
single-strandedness at the replication fork than conven- 
tional coupled leading and lagging strand DNA synthesis, 
which requires looping out of the lagging-strand template 
to permit Okazaki fragment synthesis. Alternatively, 
RNA incorporation may reflect a mechanism that minim- 
izes the deleterious consequences of collisions between 
transcription and replication complexes (30), by inhibiting 
transcription on replicating molecules (31), or restraining 
the replisome. 

The in organello mtDNA replication system provides a 
test-bed to study the role of proteins implicated in 
mtDNA replication, such as the mitochondrial RNA 
polymerase POLRMT (32), mitochondrial transcription 
factors (33), members of the mTERF family (34) and 
the DNA helicase Twinkle (C10orf2). The in organello 
system can augment in vitro systems based solely on 
purified proteins acting on artificial templates (35,36). In 
addition, for replication factors known or suspected to be 
shared between mitochondria and the nucleus, use of the 
in organello system allows physiological effects 
in mitochondria to be studied specifically, which would 
be considerably more difficult to accomplish in whole 
cells. 

The in organello labeling procedure also has the poten- 
tial to answer many outstanding questions about mtDNA 
replication in pathological states. For example, by 
studying the relative rate of incorporation of labeled nu- 
cleotide precursors into newly synthesized mtDNA, it 
should be possible to confirm or refute the idea that 
pathological Twinkle helicase mutants cause replication 
stalling (37). Nucleotide balance is also recognized as a 
critical parameter influencing mtDNA replication, as im- 
balances precipitate mtDNA depletion and disease 
(38,39). This too may lend itself to analysis in organello, 
although differences in the relative rate of uptake 
between different nucleotides will need to be taken into 
account. 



SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online: 
Supplementary Figures 1-8. 
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